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1. Introduction 



The calculation of the wave-function of a fully interacting many-body quantum system 
is a formidable challenge pQ. Hence, feasible approaches to the quantum many-body 
problem are very important if one is interested in the properties of complex multi- 
particle systems. One such approach is Nobel-Prize-winning density functional theory 
(DFT) (TJ E]. The main theorem of DFT is the Hohenberg-Kohn theorem [3j, which 
proves that the external potential of an interacting iV-body system uniquely defines the 
one-particle density of the ground-state via minimization of the corresponding energy 
functional. The mathematical foundations of DFT were extensively investigated and a 
rigorous formulation is available [U [5] . Substantial extensions of the ground-state theory 
can be found in the literature [2], e.g. to excited states or to relativistic systems. 
The mainstay of applications of the minimization principle of DFT is due to Kohn and 
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Sham [B|. The so-called Kohn-Sham scheme uses an auxiliary system of noninteracting 
particles which has the same energy and one-particle density as the interacting system. 
To make contact to the physical system one introduces the so-called exchange-correlation 
energy functional [2]. It accounts for the difference between the combined kinetic and 
interaction energy of the interacting system and the kinetic energy of the noninteract- 
ing Kohn-Sham system. The variational minimization of the energy with respect to 
the density leads to a set of coupled, nonlinear single-particle differential equations, 
the so-called Kohn-Sham equations. It can be rigorously proven that a self-consistent 
solution of these equations will generate the exact one-particle ground-state density of 
the corresponding interacting system [2] . In practice, however, the exchange-correlation 
energy functional is not known and has to be approximated. 

An exact extension of DFT to time-dependent systems was given in [7j by Runge and 
Gross. By assuming the external potentials to be Taylor expandable in time about t = to 
they could prove a one-to-one correspondence between time-dependent densities and ex- 
ternal potentials. However, the straightforward extension of the Kohn-Sham scheme to 
the time-dependent shown in [7] led to the so-called symmetry-causality paradox 

[8] . This flaw in the time-dependent extension of DFT was soon realized to be connected 
to the naive application of the usual variational principle of time-dependent quantum 
mechanics to TDDFT [9]. Only with an extension of the Runge-Gross theorem by van 
Leeuwen in \X0\ one was able to justify a time-dependent Kohn-Sham scheme. There 
it was shown that by successive solution of Sturm-Liouville boundary value problems 
one can formally construct a unique effective potential governing the time-evolution of 
the noninteracting system such that it reproduces the interacting one-particle density. 
In view of this theorem one does not need a variational approach analogously to time- 
independent DFT and can obtain the exact one-particle density via propagation of the 
time-dependent Kohn-Sham equations. 

Note that, although the extended Runge-Gross theorem [TO] shows the uniqueness of 
the effective potential, the conditions of existence were not investigated. 

The intention of this work is not to give an introduction to TDDFT, for this we re- 
fer to [XT], but to consider the mathematical foundations of the theory. In contrast to 
DFT a mathematically rigorous formulation of TDDFT is missing. Here we will take 
a first step towards this goal. We will give conditions for the existence of a solution 
to the Sturm-Liouville boundary value problems at hand. We will introduce the set of 
external potentials and time-dependent densities under consideration. We will prove 
that all orders of the Taylor expansions of the effective potentials exist, if the initial 
configurations and the different two-particle interactions as well as the external poten- 
tial of the interacting system are spatially infinitely differentiable. Given that these 
conditions hold, we will show the set of f-representable time-dependent densities being 
purely determined by the initial one-particle density and its first derivative in time at 
t = t . 
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Section 2 summarizes the basic theorems of TDDFT so far. Section 3 outlines fun- 
damental properties of the potentials and investigates the existence of general solutions 
of the corresponding Sturm-Liouville problems. In section 4 we introduce a certain set 
of potentials and f-representable densities and prove the existence of the solutions of 
all considered Sturm-Liouville boundary value problems, and hence all orders of the 
Taylor expansion of the effective potentials, under certain restrictions. Assuming these 
conditions we can evidence properties of the set of t>-representable densities. Finally we 
conclude in section 5. 



2. Time-dependent density functional theory 

We consider a general many-body Hamiltonian in atomic units of the form 

H(t)=f + V int + V([v];t). (1) 

The operator f = J2a I <Prif>l(r) (-|V 2 ) if> (r) is the kinetic term, Vi nt = 
|52ct<t' / / d 3 rd 3 r'vi nt (\r — r'\)^\(r)^\,{r')^ a i(r')'ip a (r) is the interaction and V"([V|;t) = 
/ d 3 r v(r,t) ^ a {r)^ a {r) is the external potential, where V>£(r) and tp a (r) are the 
creation and annihilation operator with spin a [12J. The interaction potential fmt(|^ — r'\) 
is arbitrary but will usually be chosen to be equal to the Coulomb interaction. The 
external one-particle potential v{r, t) typically consists of a static part, e.g. the attractive 
Coulomb potential of a fixed nucleus, and a time-dependent part, e.g. a laser pulse in 
dipole approximation and length gauge. The time-dependent one-particle density is 
defined by the expectation value of 

H r ) : = 2$(0$<t(0 ( 2 ) 

with the time-dependent density matrix p{t), i.e. 

n(r, t) = (n(r)) = tr {pit) h{r)) . (3) 
With the current- density operator 

^ := Yi^ {^( r ) V ^( r ) - [ V ^( r )] Mr)} (4) 

it is straightforward to find the usual continuity equation via application of the 
Heisenberg equation 

d t n(r,t)= -i(\h(r),H(t)]^ = -V-j(r,t), (5) 

where d t = d/dt and [•,•]_ is the usual commutator. The time-derivative of the current- 
density leads to the local force balance equation [11] 

d t ju{r,t) = -n{r,t)d u v{r,t) - d, (^(r)} - (w v {r)) (6) 

where we made use of the Einstein summation convention, i.e. summing over multiple 
indices, the momentum-stress-tensor 

TAr) := i £ | (0„#(r)) dj a {r) + (&#(r)) dj a {r) 
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-2^(^(0^(r))|, (7) 
and the divergence of the interaction-stress-tensor [13] 

WAr) : =5Z/ rfV (^int(|r-r'|)) ^(^^,(0^(04(0. (8) 

a,a f 

Now we will shortly sketch the idea underlying the Runge-Gross proof. Assume 
two external potentials v(r, t) and v'(r,t), both Taylor expandable about the initial 
time t = to, which differ by more than a merely time-dependent function c(t), i.e. 
v(r,t) — v'(r,t) 7^ c(t). If we evolve an initial configuration p(t ) = p of a finite multi- 
particle system in time with the two different external potentials, we can investigate the 
difference of the time-derivatives of the current-densities at time t = to via application 
of equation ([6]). If the time-derivatives of the current-densities differ for some order, 
then the corresponding densities will be different after an infinitesimal time-step. This 
leads to the Runge-Gross theorem [7J: 

Theorem 1 For every single-particle potential v(r,t) which can be expanded into a 
Taylor series with respect to the time coordinate around t = t , a map G : v(r,t) i— > 
n(r, t) is defined by solving the time- dependent Schrddinger equations with fixed initial 
configuration p(t ) = po and calculating the corresponding density n(r,t). This map can 
be inverted up to an additive, merely time- dependent function in the potential. 

Now let us fix the additive merely time-dependent function in the potential to be equal 
to zero by the boundary condition v (r, t) — > for |r| — > oo. Hence we have an invertible 
mapping G which obviously depends on the initial configuration po- The proof is not 
restricted to Coulombic interactions and can be applied to any reasonable interaction. 
The domain and the range of this mapping is not further investigated in the original 
paper [7J. With the mapping F : v(r, t) \— > pit) defined by the solutions of the associated 
Schrddinger equations we further find via FoG~ x : n{r, t) t— > p{t) that every expectation 
value of an operator O, i.e. O = tr[Op(i)], is uniquely determined by the density 
alone. Although this theorem holds also for noninteracting systems, it is not clear 
that every density subject to the Runge-Gross theorem in an interacting system can be 
reproduced by some effective Kohn-Sham potential in a noninteracting system. In other 
words, it is not known if the interacting v -represent able density is also noninteracting 
w-representable. Initial attempts to construct such connections in the original paper [7J 
led to the symmetry-causality paradox [8]. 

In order to overcome these problems the extended Runge-Gross theorem was introduced 
in [10]. A sketch of the proof reads as follows: 

We apply the continuity equation ([5]) to the local force balance equation ([6]), leading to 
J^n(r, t) = V ■ [n(r, t)Vv(r, t)} + (d v (d^ u {r) + W v (rj) \ . (9) 

=:q(r) 

If we now assume the external potential v (r, t) as well as the density n{r, t) to be analytic 
about t = to the different orders of the Taylor expansion are connected via equation §§§ 
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leading to [TT] 



n 



(fc+2)(r)= g (*)( r )+^ 



1=0 



k 



k 



) 



V- [n {k - l) (r)Vv {l \r)] , 



(10) 



for k > 1 where we used nP^ir) = d k n(r,t)\ t=to and q^ k \r) is defined by applying the 
Heisenberg equation /c times to g(r) at time t = t . Hence, q^ k \r) contains terms v®(r) 
up to order I = k — 1. For a second system with Hamiltonian 



and initial state p'(to) = Po anc ^ accordingly redefined operator q'(r) we can rederive 
equation ffTOj) for the primed system. Given the time-dependent density n(r,t) of the 
unprimed system we can define the potential v'(r,t) leading to the same density in the 
primed system in terms of its Taylor expansion in form of a Sturm-Liouville problem 



i=o \ / 

With this we can state the extended Runge-Gross theorem |10j : 

Theorem 2 For a Hamiltonian H{t) with an analytic potential v(r,t) about t = to we 
assume the density n(r, t) generated via propagation of the initial configuration p to be 
analytic about t — to as well. For a second system H'{t) with initial configuration po 
subject to the conditions 



finite, the analytic potential v'(r, t) leading to the same density n(r, t) is uniquely defined 
up to a purely time- dependent function. 

Again we can fix the purely time-dependent function of the potentials to be equal to 
zero by choosing the boundary condition v'(r,t) — ► for |r| — > oo. The existence of this 
potential v'(r,t) is investigated in [9] by minimization of a corresponding functional. 
However, a rigorous proof of existence was not given. 

Note, in the above schematically depicted proof of the extended Runge-Gross theorem 
one made use of the knowledge of the density n(r,t). However, we want a theory 
predicting the density. That TDDFT is a predictive theory can be seen if we use equation 
(flOl for the primed and the unprimed system and assume both to yield the same density 
n(r,t). With the definition 



H\t)=T + VU + V\[v'} ] t) 



(11) 




(12) 




v'(r, t) = v(r, t) + v a(t, t) 
we thus infer the Sturm-Liouville problem 



(14) 



V • [n (0) (r)Vvi fc) (M; r)] = (W(r) 




(15) 
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Via equation (1151) one can generate v ,<yk \r) which then can be used in equation (1101) 
to find the corresponding higher terms of the Taylor expansion of n(r,t). In return, 
the higher terms of the density Taylor expansion again determine the next term in the 
Taylor expansion of the potential. Hence, only the initial configurations and the external 
potential of the unprimed system are needed to generate the expansion of the density 
about t = t . However, the special form of g"(r) defined by successive application 
of the Heisenberg equation led to the question wether a simultaneous solution of the 
unprimed system is required in principle [HI [15] . 



3. Properties of the potentials and the Sturm-Liouville problem 

In order to define the sets of one-particle potentials and the associated one-particle 
densities, i.e. the v -represent able densities, we start by looking at the properties of 
the corresponding Hamiltonian H(t). We demand the Hamiltonian to be self-adjoint for 
every time t on the domain of the kinetic energy operator dom(T). The free Hamiltonian 
Hq — T + Vint with Vi n t the Coulomb interaction can be shown to fulfill this constraint 
[5] and any other interaction V( nt under consideration will be assumed to do so as well. 
This condition is trivially fulfilled if we choose V{ nt = 0. With the theory of Kato 
perturbations [5] we find the simple condition 

v(t) G L 2 (M 3 ) + L°°(IR 3 ) (16) 

for every time t. Here L 2 and L°° are the usual Lebesgue quotient spaces with norm 
||.||2 and ||.||oo; respectively. L 2 (IR 3 ) + iL°°(M. 3 ) is a Banach space with the norm 



inf{|K(t)|| 2 + |Mt)||oo \v x {t) G L 2 (R 3 ), v 2 (t) G L°°(M 3 ), 

v{t) = v 1 {t)+v 2 {t)\. (17) 



The effective potentials v'(r,t) leading to the same density as in the Coulombic 
system exists in accordance to the extended Runge-Gross proof if all orders of equation 
( fl2l) have an existing solution. Accordingly, this is true if all orders of equation ( {15}) 
have a solution with v'(r,t) = v(r,t) + v&{r,t). A way to investigate the existence of 
such solutions to such Sturm-Liouville boundary value problems is within the following 
framework. 



Let us first introduce the bounded open domain Q C M 3 with piecewise C 1 boundary 
dfl. For the real Hilbert space L 2 (Q) with scalar product denoted as (•, -) 2 and norm 
|| ■ || 2 we define the Sobolev space [5] 



W 1)2 {VL) = jw G L 2 {Vl) d jU G L 2 {Q), j = 1,...,3} (18) 

with norm ||«||i 2 = IMIi + ll^^lll- Here the partial derivatives djU are understood 
in the distributional sense. Further the real Hilbert space Hq(Q) is defined to be the 
closure in W 1)2 (Vt) of the infinitely differentiable functions compactly supported in Q. 
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Hence, functions u G Hq(Q) fulfill the boundary condition u = on <9f2 naturally. The 
previously considered Sturm-Liouville problem fll5D in shorthand notation reads as 



V ■ (nVv) = C (19) 

and defines a bilinear form Q on Hq(Q) by 

Q(u, v) = (u, -V • (nVv)) 2 = (Vu, nVt>) 2 (20) 

where we used integration by parts and that functions in Hq(Q) vanish at the border. 
We take advantage of the fact that ffl9|) has a (weak) solution v G Hq(Q) if, and only 
if, Q(u,v) = —(11,02 f° r & U M ^ -^o(^)- This immediately leads us to the necessity 
( G L 2 (Q). Now the answer to the question of solvability is at hand with the Theorem 
of Lax-Milgram. [16] 

Theorem 3 (Lax-Milgram) Let Q be a coercive continuous bilinear form on a Hilbert 
space TC. Then for every continuous linear functional f on 7i, there exists a unique 
Uf G TC such that 

Q(u,u f )=f(u) (21) 
holds for all u ETC. 

A bilinear form Q is said to be coercive if there exists a constant c > such that 
Q(u,u) > c||w|| 2 for all u G 7i. In our case this can be established by means of the 
Poincare inequality 

IMh < A||Vu||a, Vw G H%(0) (22) 

where < A = X(Q) < 00. As an additional assumption we add that n is bounded by a 
constant m > almost everywhere on Q from below. Then 

Q(u, u) = (Vu, nVu)2 > m\\ Vu||| (23a) 
X 2 Q(u,u) > A 2 m||VM||2 > m||u||| (236) 

Combination of these results yields 

Q(u,u) > + ||Vu||2) = ^LjUull^. (24) 

We thus have established the coercivity of Q. As for the continuity we add another 
assumption on n that is boundedness from above by a constant M > almost 
everywhere on Q. 

\Q(u,v) - Q(u ,v)\ = \Q(u - u ,v)\ = \(V(u - u ),nVu) 2 \ 

< M||V(w - uo)|| 2 ■ ||Vv|| 2 < M\\u - w ||i )2 ■ |M|i,2 < 00 (25) 

These restrictions on n also imply that the differential operator defined by the left 
hand side of ( fT9l) is elliptic. If n can be assumed to be continuous on the closed domain 
fl, then n also attains its extremal values on Q and the restrictions reduce to the form 

< n < 00 on H. (26) 
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Now everything left to show for the application of Theorem [3] is that the right hand 
side — (-,02 is indeed a continuous linear functional on the real Hilbert space Hq(Q). 
This is easily established by considering for arbitrary u, uq G Hq (O) 

|(m,C) 2 - (wo,02| = \(u-UoX)2\ 

< \\u - uoh ■ HClla < \\u - II 1,2 • HClla < °o. (27) 
We subsume our results for the solvability of the Sturm-Liouville problem ( fl9i) in 
the following corollary. 

Corollary 1 Consider the Sturm-Liouville problem V ■ (riVv ) = ( on the bounded open 
domain Q C M 3 with piecewise C 1 boundary. Let £ G L 2 (f2) and n : f2 — > R 6e almost 
everywhere bounded byO<m<n<M<oo. Then there exists a unique solution 
v G /^(ft). 



Looking back to the original formulation of the Sturm-Liouville problem ( 1151) we 
note that if n^ ' fulfills the prerequisites of n in Corollary [1] and the right hand side 
is indeed in L 2 (Q) then all orders of the potential are uniquely defined in Hq(Q). 
We now want to examine which constraints on the higher orders are sufficient for 

e L 2 (Q). For this we assume q^ and q'^ already in L 2 (f2) for all k and thus we get 
a unique G #q(^) trivially by Corollary [U For k > 1 we apply inductive reasoning: 
Let us assume G Hq(Q) be given uniquely for / < k. Then ffl5l) in a shorter notation 
reads as 



V • (n^Vvg ) = ff « - - ^ ( k j V • (n^VcJ) . (2 



The (elliptic) differential operator on the left defines the same coercive continuous 
bilinear form on (Q) as before, thus everything to show for an application of Theorem 
[3] is that the right hand side yields a continuous linear functional. The sum of such 
functionals is again linear continuous therefore we can examine all terms separately. 
For g( fc ) and q'^ the same reasoning as in (j2"Tj) applies. Finally with u, u G Hq (fl) 
continuity of the separate terms of the sum is established by 

(u, V ■ (>-<) V^) } 2 - (u Q , V ■ (nM V «J) } 2 

< \\u - Uo||i,2 ■ Wn^Vv^h < oo if n {k - l) Vv { 2 G L 2 (fi). (29) 

By induction we know that G Hq(Q) and thus Vf^ G L 2 (f2) so bounded 
almost everywhere by some > for all A; is a sufficient condition for Theorem [3] to 
be applied. We subsume this in a second corollary. 

Corollary 2 Consider the system of Sturm-Liouville problems l[T5\) on the bounded open 
domain f2 C R 3 with piecewise C 1 boundary. Let gW,<p fc ) £ L 2 (f2) and n( fc ) : — > M 
fre almost everywhere bounded by nS k ' < with > for all k and additionally 
n (o) > m > almost everywhere. Then there exists a unique sequence of solutions 
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If we turn to the problem of a classical solution we refer the reader to the Weyl 
Lemma as given in [T71 IE] in various forms. (In a more general formulation it can be 
found as the "fundamental theorem on weak solutions" in [19J.) There the operator 
Kv = V • (nVv) has domain 

dom(K) = |m u G C\(}i), u G C 2 (n), Ku G L 2 (n); w = on <9fij . (30) 

Theorem 4 For n > on Q and n G C 3 (f2) and the boundary condition v — on dQ, 
the equation 

V • (nVo) = C (31) 

has a classical solution v G dom(iT) if ( G C 1 (fi) or C is Holder continuous, i.e. 
ICO") — C( r ')l — ^l r — r '\ a f or a M r i r ' ^ ^ with h and < a < 1 independent of 
r, r'. 

Wether C (fe) 

in our actual problem (fT5l) fulfills one of the conditions for a weak or 
classical solution, respectively, depends on the properties of the initial configurations 
and on the interactions under consideration. The terms and q'^ implicate already 
for k — spatial partial derivatives of order 4 and spatial partial derivatives of the 
involved interaction potential of order 3. Hence, to have a well defined Sturm-Liouville 
problem, the wave-functions of the initial configurations and the interaction potentials 
have to fulfill certain restrictions with respect to their spatial behavior. 

4. Sets of potentials and w-representable densities 

We will introduce the set of external potentials for the extended Runge-Gross theorem 
in accordance to the classical Sturm-Liouville theory. Therefore, the defined sets will 
only be subsets of the actual sets of v-representable densities and potentials connected 
via the extended Runge-Gross theorem. Further, we will restrict our considerations on 
the above introduced domain Q C M 3 . We assume the boundary to be far away from 
the center of the system such that it will not influence the dynamics. One has to be 
careful at this point as we will assume the initial one-particle density to be nonzero on 
Q, hence also on the boundary dQ. This is different to the usual notion of a physical 
system restricted to a finite region, where one assumes an infinite boundary potential 
to restrict the wave-function to this domain. In this case the wave-function and thus 
the one-particle potential will be zero at the boundary. 

For a free Hamiltonian Hq = T + Vint assumed self-adjoint and an initial configuration 
p at time t = t we have 

V(/3o, Vint) := v analytic about t = to,v(t) G dom(i^), 

v real; n[v] analytic about t = t for p and Vi nt |. (32) 

Here n[v] is the time-dependent density, defined via the propagation of the initial 
configuration p with the Hamiltonian H(t) = T + Vi nt + ^([u]; t). It is straightforward 
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to proof self-adjointness of this Hamiltonian by application of the Kato perturbation 
theory [5] as one can use v(t) G L°°(Q). Further we define the set of f-representable 
variations by 

SAf(poAnt) := {Sn 5n(r,t) = ^ ^ k \[v}; r)(t - t ) k for p and V int , 

k=2 

ve V(p ,Mnt)}. (33) 

The set of v-representable densities is an affine set 

AT(p ,^nt) := n^(r) + n^(r,t) + 6Af(p ,V int ), (34) 

where n^°'(r) = n(r,to) is the initial density and n^(r,t) = tr(poV ■ j(r))(t — to) i n 
accordance to (113 ah and ( 11361) . respectively. For these sets we then have in accordance 
to the Runge- Gross theorem an invertible mapping 

Vf> :Af(p ,V in t)^V(p ,V int ) (35) 
n(r,t) ^Vp ([n];r,t) 

connecting the w-representable one-particle densities with the external potentials. Never- 
theless, if we now define a second mapping for a different initial configuration p' subject 
to the conditions (113al) and (113 61) . and a different interaction V{ nt , we do not know if n is 
simultaneously element in M(po, Vint) and Af(p' , V^ nt ). This, however, is of fundamental 
importance if we want a rigorous formulation of the time-dependent Kohn-Sham scheme. 

To achieve this goal we introduce further restrictions. We will assume smooth inter- 
actions and initial configurations in what follows. This excludes the usual Coulombic 
interaction as it is not infinitely differentiable at the origin. However, one may regu- 
larize the Coulombic interaction by a so-called soft-core interaction, i.e. by replacing 
\r\ — > y/r 2 + e and e > 0. 

Due to equation (fit)]) we have a direct connection between f-representable densities and 
potentials. Hence, we can formulate the following lemma 

Lemma 1 Let p be chosen such that all its wavefunctions are in C°°(Q N ), 
v Taylor expandable about t = t , v^ k ' G C°°(n) V k, and v- mt (\r — r'\) infinitely 
differentiable. Then n^ k \r) G C 00 (n) for all k. 

Proof. We will use equation (flOl) . Obviously we have n^(r) and n^(r) in C°°(f2). 
Thus n^ 2 \r) is in C°°(f2) if q^(r) is infinitely differentiable, where 

g(°)(r) = tr 



(36) 



q(r) consists of partial derivatives with respect to r and of derivatives of tw(l r — r '|)- We 
have assumed v- mt (\r — r'\) infinitely differentiable. Hence, we have q(°'(r) G C°°(Q,). For 
72,( 3 )(r) we need to know q^'(r). This is the commutator of q(r) with H(t) at t = to. All 
functions in H(t) are infinitely differentiable. Again the above reasoning applies, and 
we find q( l >(r) G C°°(f2). All higher terms are to be found via successive application of 
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the Heisenberg equation for q(r) with H(t) at t = t . The only difference to the above 
reasoning is the appearance of v^ k ' (r)-terms, which are again infinitely differentiable. 
Therefore one can successively construct all rS k \r) G C°°(fi). □ 

Now we introduce the restricted set of smooth one-particle potentials 

V*(f>o,V int ) = {v\v eV(p ,V int ), v(t)eC^(Q)Y (37) 
and by lemma [1] the corresponding smooth f-representable one-particle densities 

oo _^ 

N*(po,V int ) = jn n(r,t) = 7J n (k \[v}; r)(t - t ) k for p and V int , 

k=0 



veV*(po,Vte)\. (38) 

With this we can reformulate the extended Runge-Gross theorem as follows. 

Theorem 5 Let p and VJ nt be infinitely differentiable, n(r,t) G A/"*(p , Vmt) an d 
Vp ([n};r,t) = v(r,t) G V*(p ,Vi n t) the associated external potential. For a system 
with infinitely differentiable interaction V( nt and the initial configuration p' Q consisting of 
infinitely differentiable functions subject to the constraint 

n(r, t ) = n (0) (r) = n'{r, t ) > 0, (39) 
tr(p V-j») = tr(p' V-i(r)), (40) 
there exists a unique effective potential depending on both initial configurations 

oo , 



k\ 

k=0 



where v^\r) is defined via 



n (0) (r)W^(r 



k-i 



g (fc)( r ) _ g '(*)( r ) _ J2 ( * j V • [n^ k ' l \r)Vv' W (r)] , (42) 

with v ' = (y + Up ,^) G V*(p' Q , V{ nt ) generating the same density. It holds that 

A/-*(p ,VW =M*{p',XJ =Af*(n(°\nV). (43) 

Proof. From the proof of lemma [T] we know that all q^ k ' and cp" are infinitely 
differentiable. As we have assumed no(r) > we can apply theorem H] from which 
it is clear that 

V- [n (0) (r)W (fc) (r)j = 



q {k \r)-q'^(r) - 



k-l i , 

K 



J2[ j V- [n' {k - l \r)Vv'^{r)] (44) 



1=0 



has an existing solution for k = if the right hand side is C 1 (f2). Obviously iA°)(r) 
exists due to theorem H] and is infinitely differentiable. In the next step we can use 
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iA°)(r) in the Sturm-Liouville equation denning i>'W(r). Again existence is guaranteed 
and we have u'W(r) G C°°(f2) [20]. One can now successively construct Vp 0t p . Then 
{v + Vp (h p' a ) is given via its Taylor series within its radius of convergence in accordance to 
the extended Runge-Gross proof [10]. This construction holds for every n G J\f*(p , V int ), 
and we have n e J\f*(p'o, ^int) as we ^- Hence, the set of w-representable densities does 
not depend on the smooth interaction or on the smooth initial configuration. □ 

Here it became obvious why we restricted our considerations to infinitely differentiable 
initial configurations and potentials. With this assumptions we can guarantee the exis- 
tence of all the classical Sturm-Liouville boundary value problems on Q. For the general, 
i.e., weak case, we need to make sure that all q( k \q'( k *> and hence are in L 2 (f2), in 
order to proof the existence of a solution using corollary [2j 

The special case of a rigorous Kohn-Sham theorem is straightforward as V- mt = is of 
course infinitely differentiable. One finds that for the above restrictions all interacting- 
v -represent able densities are noninteracting-t>-representable because A/"*(p , Vint) = 
jV*(po, 0). Only a noninteracting initial configuration is needed. The condition of 
7T,(°)(r) > for the existence of the effective potential may be relaxed if there exists 
some time t\ in a sufficiently small neighbourhood of to for which n(r, ti) > 0. Then we 
could use p(t\) as new initial state and prove existence at that time provided we also 
have the corresponding p'{t\). 

5. Conclusion 

Under certain assumptions we can state sufficient constraints on the one-particle 
density such that the existence of the effective potential, possibly in the weak sense, is 
guaranteed. However, only for classical solutions of the corresponding Sturm-Liouville 
boundary value problems we can reformulate the extended Runge-Gross theorem such 
that existence of the effective potentials is granted. As long as we consider smooth 
initial states and smooth interactions the Kohn-Sham system exactly reproduces the 
physical one-particle density. In general, as pointed out in [9], it seems safe to assume 
existence of the Kohn-Sham potential for physical systems. Nevertheless, a rigorous 
proof of principle is of importance for the foundations of the theory. 
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